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The 3X + 1 Problem 


The 3X + 1 problem was devised in 1967 by Prof. 
Richard Andree of the University of Oklahoma, to demonstrate 
the concepts of recursion and algorithms to a class of 
high school students. The problem is this: for a 
positive integer, N, let X equal N and proceed as follows: 


Replace X by X/2 if X is even. 
@ Heplace X by 3xX+1 if X is odd. 
Stop when X equals 1. 
Call the number of terms so generated A, 
counting the original number. 


Thus, for N = 9, the pruceated sequence is 9, 28, 14, 7, 
22, Ml, Bh tye Ges ea, sh lon can len So Ne, oy, 4 25 a. 
and A = 20. 


Notice that in performing this calculation, we have 
certified the convergence of many other values of N, and 
also determined their A values. 


The first natural question that arises is, will 
the process converge for all values of N? This seems to 
be difficult to prove. Convergence has been established 
by computation for all values of N up to 200,000,000, 
utilizing the following shortcuts: 


1. Only odd values of N need be tested, if the 
computation proceeds upwards in the N values, since for 
any even N the half value will then already have been 
tested. 


2. It is not necessary to proceed to X = 1 to prove 
& convergence, but only to a value of X less than N. 


3. Numbers of the form 4K+1 will follow this 
sequence: 


4K+1 
12K+4 
6K+2 
3K+1 


and hence will always converge, provided that values of N 
are tested in increasing order. Thus, only odd values 
of N of the form 4K+3 need be tested. ‘The table of 
Figure 1 shows the slow growth of values of A. 


The A values seem to cluster in a peculiar way. 
For example, these values of N all have an A of 39: 


105,614, 612, 613, 614, 
628, 629, 630, 631, 632. 


For any limited range of N, there will be comparatively 
few values of A. Thus, in the range of N from 90,000 
through 94,999, over half of the A values are these: 
SE, Ul, He, Dee, eis) 194s). 


An N can be found for any given A, JES), it eye 


aay a fa 
AS Teo WEBS NE ak ys lft values! of N are restricted to 
being odd, then not every A can be produced, and empirical 
studies indicate that relatively few A values exist. 
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23529 
26623 
34239 
35655 
52527 
77031 
106239 
142587 
156159 
216367 
230631 
410011 
Dds 
626332 
637799 
1117065 
ANCHE Sy 
L/AS ILS 
2298025 
3064033 
3542887 
37 32423 
5649499 
6649279 
8400511 
11200681 
14934241 
15733191 


The table shows 
Only odd values 


the first appearance of successively larger values 


of A, with the corresponding N. 


Results from the 3X + 1 problem. 
of N were being examined. 


Figure 1. 
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What remains to be done with the 3X+1 problem? 
Values of N greater than 200,000,000 could be tested 
for convergence (although this path is losing its charm). 
One new computer approach would be the following. Clear 
to zero a block of bits each of which represents an odd 
integer. Now, in testing an odd N for convergence, 
store a 1 at the bit position corresponding to any odd 
number greater than N that is met along the way. Thus, 
in testing N = 9, the bit positions corresponding to 
values of 11, 17, and 13 would be set to l. The 
remaining zeros in the block of bits would then point 
to the next value of N that needs to be tested. The 
shortcuts listed earlier can still be applied. 


It might be concluded that the process will probably 
always converge, and the proper avenue of attack on the 
convergence problem will be analytic. As a computer 
problem, work needs to be done on the Gistribution of 
the A values. 


woddd 
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fhe Wells/Ulam Conjecture 


Writing in the September, 196/) Sclentific American, 
Stanislaw Ulam reported on a conjecture made by Mark 
Wells, to the effect that multiples of 3, expressed in 
binary notation, should most often have an even number 
of l-bits. This is stated as a theorem: "Among all 
the integers divisible by three from 1 to 2N, those that 
have an even number of 1's always predominate, and the 
difference between their number and the number of those 
with an odd number of 1's can be computed exactly: it 
is 3(N-1)/2 0 


The table shows actual counts up to 216 | plus two 
more observations that are not at powers of 2. The 
value given by the formula is rounded to the nearest 
integer. There is also shown the ratio of odds to 
evens, which seems to tend toward 1.000 as N increases. 
The repetition in the ratio indicated by the brackets 
on the far right is probably significant. 
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65536 13109 8736 
500000 92463 #74203 
1000000 184900 148433 
1500000 270475 229525 


Ratio of odds to evens 


> {Nei )/2 
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Desk Calculator Review 


Hewlett-Packard HP-35 


The HP-35 is unquestionably the leader in the pocket 
calculator field; it is the machine against which all 
others must be conpared. It is in a class by itself, 


both on price ($400), its precision (10 significant digita), 


and its functions (natural and common logarithms, 
exponential, power function, direct and inverse sine, 
cosine, and tangent). A 4-level stack storage allows 
elaborate cascading of operations and forces the use of 
Polish (Lukasiewicz) notation. Thus, for the operation 
A/B = Q, most calculators follow the sequence: enter A, 
press divide, enter B, press equalia. On the HP-35, 

the sequence is: enter A, enter B, press divide. The 
machine also has, in addition to the 4-level stack, 

one additional storege word. 


Generally, the claim for 10 significant digit 
precision holds good. Inverse operations such as 


X, reciprocal, resiprocal = XK 


ie = (Gls in. @ = Xx 
return to the exact original vaiue for a wide range on X. 


All arithmetic is floating point in the range from 
107° to 1040 and outeide that range is in scientific 
notation, up to 9.999999999E99. Illegal operations, 
such as division by zero or the logarithm of a negative 
number cause the display to blink. The sequence of 
operations: enter 2, square root, square root, square 
root, square root, square root, square, square, square, 
square, square, equals 2.000000022 indicates the inherent 
accuracy of the algorithms that have been programmed 
into the machine. On the other hand, the sine function 
appears to apply an algorithm directly without reduction 
by multiples of two pi; thus, sine 30 = .5 but sine 750 = 
-4999999986. Similarly, sine 720 = 4k-09, 


Early models of the HP-35 had an error in one of 
the LSI chips. in a mailing to registered owners, the 
company stated "In most cases, the error range in the 
answers is from one one-hundreth (sic) of one percent 
to a maximum of one percent." The first such error 
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listed in that mailing is: 


5.729577893E-3 (HP~-35) 
01145916 (true) 


arcsin .0002 
arcsin .0002 
which is a 50% error, 


In any event, this set of bugs has been corrected on 
later serial numbers, and Hewlett-Packard has agreed to 
modify the early machines. 


The HP-35 is a BEST BUY by any criterion. No 
user, after a few weeks of use, would be likely to want 
to go back to a simpler machine at any price. The 
machine represents a significant jump in calculator 
capability and, after a year of sales, has no competitor. 


THE CALENDAR 


The present calendar in use throughout most of the 
world began in 1582; it was adopted by Great Britain and 
the Colonies in 1752. The calendar repeats precisely 
every 400 years. The tables below show the statistics 
for any period of 400 successive years. 


Table 1 shows the distribution of month types 
for the 4800 months. There are 44 months of 28 days 
that begin on Sunday; 399 months of 31 days that begin 
on Saturday; and so on. 


Table 2 shows the distribution of days of the 
month. For example, the 13th of the month falls on 
Friday 688 times, which is more often than any other day 
of the week. 


Table 3 shows the distribution of the 14 year 
types. Thug, 44 years are not leap years and begin on 
a Thursday. Of the 97 leap years, 15 will begin on a 
Friday. 


Table 4 shows the years in the current century 
that are of each of the 14 possible types. 
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Table l: 


28 day month 
29 day month 
30 day month 
31 day month 


Table 2: 


Table 3: 


Not leap year 
Leap year 
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Table 4: 


1905 
1995 


1900 
1990 


19021 
1991 


1902 
1997 


1903 
1998 


1909 
1999 


1910 
2005 


1928 


1912 


1924 


1908 


1920 


1911 
2006 


1906 
2001 


1907 
2002 


1913 
2003 


1914 
2009 


1915 
2010 


1921 


1956 


1940 


1952 


1936 


1948 


Normal years starting 
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Normal years starting 


1917 1923 1934 1945 


Normal years starting 


1918 1929 1935 1946 


Normal years starting 


1919 1930 1941 1947 


Normal years starting 


Poe 5 t9SL e1gke 1953 


Normal years starting 


1926 1937 1943 1954 


Norwal years starting 


ese ar geO? 955 


Leap years starting 
1984 2012 

Leap years starting 
1968 1996 2024 

Leap years starting 
1980 2008 

Leap years starting 
1964 1992 2020 

Leap years starting 
1976 2004 


with Sunday 


1961 1967 1978 


with Monday 
1951 1962 1973 


with Tuesday 


1957 1963 1974 


with Wednesday 


1958 1969 1975 


with Thursday 


1959 1970) 162% 


with Friday 


1965 1971 1982 


with Saturday 


1966 1977 1983 


with Sunday 


with Monday 


with Tuesday 


with Wednesday 


with Thursday 
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1989 


L979 


1985 


1986 


1987 


1993 


1994 
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Leap years starting with Friday 
1904 1932 1960 1988 2016 

Leap years starting with Saturday 
1916 1944 1972 2000 


s 


Book Review 


Program Test Methods, edited by William C. Hetzel, & 
Prentice-Hall, 1973, 311 pages plus a bibliography of 
375 citations plus an Index of Concepts. 


This is the Proceedings of a conference on Computer 
Program Test Methods, sponsored by the University of 
North Carolina and the ACM Special Interest Group on 
Programming Languages, and held at Chapel Hill June 21-23, 
1972. it is the first attempt to bring together what 
is known about software validation. 


Program testing--how to certify a program as 
performing the way it should--ought to be a fundamental 
concept in the learning of computing. As a basic 
concept, it should be an important topic in introductory 
texts. But an examination of any dozen textbooks taken 
at random will reveal how carefully the subject has been 
avoided. The majority of texts ignore the subject 
completely. Those that do mention it usually get it 
nicely confused with debugging and then compound the 
confusion in the student's mind by offering him fatuous 
or worthless guidelines for devising test procedures. 

It would not be farfetched to say that beginning students 
are systematically being taught how to produce garbage @ 
with computers. 


To quote from the book's Preface: 
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The growth of a testing discipline has been 
painfully slow, despite the acknowledged 

need for better quality software. Even the 
acope of what is or is not a testing activity 
is not well defined. Terminology in the 
field is unclear and literature that would 
establish some foundations has not been 
written. Testing as an activity needs and 
requires more attention. 


A sense of increasing urgency seems prevalent 
and much research is underway. In general, 
however, the methodology, techniques and 
theory of program testing are entirely 
inadequate. It is hoped that this book is 
a start toward the solutions needed. 


The 23 pepers in the book are a beginning, at 
least. Each of the eight sections begins with some 
pertinent quotations, one of which is worth repeating: 


In the space of one hundred and seventy—six 
years the Lower Mississippi has shortened 
itself two hundred and forty-two miles. 

That is an average of a trifle over one mile 
and a third per year. Therefore, any cain 
person, who is not blind or idiotic, can see 
that in the old Oolitic Silurian Period, 
just a million years ago next November, the 
Lower Mississippi River was upward of one 
miliion three hundred thousand miles long, 
and stuck out over the Gulf of Mexico like a 
fishing-rod. And by the same token any 
person can see that seven hundred and forty- 
two years from now the Lower Mississippi 
will be only a mile and three quarters long, 
and Cairo and New Orleans will heve joined 
their streets together, and be plodding 
comfortably along under a single mayor and 

a mutual board of aldermen. There is 
something fascinating about science. One 
gets such wholesale returns of conjecture 
out of such a trifling investiment of fact. 


--Mark Twain (Life on the Mississippi) 


Program Test Methods is certain to become one of 
the milestones in software literature. The various 
authors addressed themselves to the stated problem; 
they suggested feasible criteria for software certification; 
a methodology to apply to future software projects; and 
the boundaries of current capability in this most 
sensitive field. 


--JdJ 
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SUBFACTORIALS 


Subfactorials are defined by the recursion: 
N 
'N = N- !(N-1) + (-1) 


with subfactorial 1 defined as zero. In the following 
table, all values up to !28 are exact; for those that 
follow, the first 30 significant digits are given, 
rounded from the 3lst digit, followed by the number of 
digits to the decimal point. 


it 0 31 302501328894 190910970370027530 004 
2 it 32 968004252461410915105184088096 005 
B 2 33 319441403312265601984710749072 007 
4 9 34 108610077 126170304674801654684 009 
5 uy 35 380135269941596066361805791395 010 
6 265 ho 300158458444475693321518926221 018 
7 1854 50 111887 196107824805046 302580708 035 
8 14833 60 30611200890 300759 3241410795996 052 
9 133496 70 44.06670251980594417 26775008254 068 
10 1334961 80 263289 31863123072598092 3677587 089 


11 14684570 90 54656435875304:08685078 36153046 108 
12 176214841 100 3433279598416 38047651959775268 128 
13 2290792932 

14 32071101049 

15 4810665157 34 


16 T697064251745 

17 130850092279654. 

18 2355301661033953 
19 447507 31559645106 
20 895014631192902121 


21 18795307 255050944540 

22 413496759611120779881 

23 951042547 1055777937262 

24 228250211305338670494289 
25 57062552826 33466762357 224 


26 148 362637 348470135821287825 
27 400579120840869 3667174771274 
28 11216215383544 342268089 359567 3 


29 3252702461227 859257 74591427452 001 
30 9758107 38368357777 323774282355 002 


